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17 SUMMARY 


18 Mesozoic was the age of reptiles that not only occupied the land and sky but also 
19 adapted to the sea. Sauropterygia, from  turtle-like placodonts, lizard-like 
20  pachypleurosaurs, and predatory nothosaurs to long-necked pistosaurs including 
21 _ plesiosaurs, is one of the most predominant lineages among secondarily aquatic reptiles. 
22 However, we know little about the early evolution of sauropterygians involving the 
23. morphological gaps among aforementioned sub-groups and their ancestral body-plan, 
24 because the Early Triassic fossil remains are scarce shortly after their origin. Here we 
25 report a skeleton from the Early Triassic in South China, representing the oldest known 
26 complete specimen’ related to Sauropterygia. It can be referred 
27. to Hanosaurus hupehensis and shows a mosaic morphology combining the characters 
28 of multiple sauropterygian sub-lineages. We constructed an updated character matrix 
29 obtaining the hitherto most comprehensive phylogeny of Triassic sauropterygians, 
30 and Hanosaurus hupehensis is stably resolved as the basal-most member of the 


31 Sauropterygiformes (clade nov.). Sauropterygians were previously considered using 
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limb propulsion for the movement that is distinguishable from ichthyosaurs and other 
marine reptiles, but this skeleton reflects an ancestral body plan of sauropterygians 
unexpectedly developing an elongate trunk and four short limbs, being different from 
most sauropterygians but more convergently similar to the basal members of other 
marine reptilian lineages in body shape. After this convergence, we confirm the rapid 
adaptive radiation of sauropterygiform reptiles following the end-Permian mass 
extinction. Our results provide an evolutionary framework of sauropterygian marine 


reptiles and highlight the integrations of both convergences and divergences when an 
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emerging animal lineage arises. 
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13 


14 RESULTS 


15 Marine reptile is a great model for evolutionary study by exhibiting secondarily aquatic 
16 adaptations!. Sauropterygia” is one of the most diversified and dominant marine 
17 reptilian groups having been studied for about two hundred years. Based on abundant 
18 complete skeletons of sauropterygians and other significant marine reptiles*+ 
19 discovered during recent two decades, especially from southern China®, we could have 
20 better glimpses of sauropterygian origin and early evolution. Here we report a nearly 
21 complete skeleton referred to Hanosaurus hupehensis®. In this research, we aim to 
22 resolve the relationship within sauropterygians with related marine reptiles such as 
23 saurosphargids, illuminate the very early history, and evaluate the adaptive radiation of 


24  sauropterygian and allied reptiles. 
25 


26 Systematic paleontology. 


oF. Diapsida Osborn, 1903 
28 Sauropterygiformes clade nov. 
29 Hanosaurus hupehensis Young, 1972 


30 Holotype. Institute of Vertebrate Paleontology and Paleoanthropology (IVPP) V 3231, 
31 a skull (Supplementary Fig. 1) in dorsal view and a partial postcranial skeleton 


32 (Supplementary Fig. 2) in ventral view’. 
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1 Type Locality and horizon. Songshugou, Xunjian Village, Nanzhang County, Hubei 
2 Province, China; dolomitic limestone (Member II) of the Jialingjiang Formation, 


Olenekian, Early Triassic°. 


4 Referred specimens. IVPP V 15911, a nearly complete skeleton in ventral view (Fig. 1, 
5 Supplementary Fig. | and 2). 


6 Locality and horizon. Yingzishan, Yuanan County, Hubei Province, China; dolomitic 


limestone (Member II) of the Jialingjiang Formation, Olenekian, Early Triassic. 


Amended Diagnosis. A medium-sized sauropterygiform reptile with the exclusive 

9 combination of snout without constriction; supratemporal fenestrae smaller than orbits; 
10 nasals longer than frontals; mandibular articulation far posterior to the level of occiput, 
11. ~with retroarticular processes well-developed; anterior teeth short and conical; posterior 
12 teeth constricted at the base, with the lingual surface of tooth crown slightly concave; 
13 cervical ribs with evident anterior processes; dorsal ribs pachyostotic; coracoid not 
14 waisted; pubis extremely round; ischium kidney-shaped; no distinct thyroid fenestra 


15 between pubis and ischium; and minimally four carpals and three tarsals. 
16 Description. 


17 This new specimen, IVPP V 15911, is entirely articulated and exposed in ventral view, 
18 _- with a total length of 79.4 cm (Fig. 1 A, B) being comparable in size to [VPP V 3231, 
19 the holotype of Hanosaurus hupehensis. The cranial shape is of an elongate triangle 
20 ~=(Fig. 1 C, D). The snout is neither pointed nor blunt without lateral constriction. A pair 
21 ~~‘ of large pterygoids meet along the midline suture anteriorly and posteriorly leaving a 
22 narrow interpterygoid vacuity. In the holotype, the preorbital and postorbital regions of 
23 the skull are nearly equal in length, and the supratemporal fenestra is smaller than the 
24 orbit (Supplementary Fig. 1), which is similar to saurosphargids®, pachypleurosaur’, 
25 and pachypleurosaur-like reptiles’ but on the contrary to derived eosauropterygians 
26 such as nothosaurs and pistosaurs*. The same positions and proportions of orbits and 
27 supratemporal fenestrae can be tentatively identified in our new specimen on the right 
28 side of the skull from the ventral view. The short mandibular symphysis is mainly 
29 contributed by the dentary. The retroarticular process is well developed as in the 


30 holotype, with a blunt termination. 


31 + ‘There are at least 6 dentary teeth exposed on the anterior portion of the right dentary 

32 and more than 10 maxillary teeth visible on the right side (Fig. 1 C, D). There is neither 

33 obvious diastema between premaxillary and maxillary teeth nor prominently enlarged 

34 fangs. All the teeth are short, sharp, and likely thecodont. The anterior teeth are 

35 generally coniform, slightly procumbent; they are relatively larger than the posterior 

36 teeth and possess longitudinal striations. The lateral teeth are smaller, mildly overbite, 
3 


bioRxiv preprint doi: https://doi.org/10.1101/2022.04.25.489368; this version posted April 26, 2022. The copyright holder for this preprint 
(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is 
made available under aCC-BY-NC-ND 4.0 International license. 


and show relatively smooth surfaces of the crowns with slight constrictions at the bases 
of the crowns (Fig. 1 C, D; Supplementary Fig. 2). Additionally, same as those in the 


1 
2 
3 type specimen (Supplementary Fig. 1 and 2), the crowns of the lateral teeth are slightly 
4  lingually concave showing leaf-shaped outlines, which resemble those in 
fs) 


saurosphargids®!°. 


There are 15 cervical, 25 dorsal, and 3 sacral vertebrae, and at least 16 caudal vertebrae 
in IVPP V 15911. The number of cervical vertebrae here in Hanosaurus (15 in total 
with 13 anterior to the pectoral girdle) is fewer than all Triassic eosauropterygians. We 
identify the last cervical vertebra according to the corresponding rib not articulated with 
10 the sternum!'! instead of the position anterior to the pectoral girdle. Therefore, the 
11 previously reported cervical counts’? should be revised as about 20 in Neusticosaurus 
12 and 18 in Serpianosaurus based on our first-hand observation, which develops fewer 
13 cervical vertebrae than other Triassic eosauropterygians but more than Hanosaurus. 
14 Including the ribs, the entire vertebral column is conspicuously pachyostotic. The 
15 centrum is wider than long and anteriorly narrow and posteriorly wide in ventral view 
16 ~=©(Fig. 1 A, B, G, H). The cervical rib possesses a distinct anterior process in vertebrae 2 
17 ‘to 13 (Fig. 1 E-H; Supplementary Fig. 2 I-L). The last cervical rib is elongate, single- 
18 headed, and nearly as long as the first dorsal rib, but its shaft is thinner than the latter 
19 (Fig. 1 G, H). The anterior caudal ribs are anteroposteriorly broader than the gap width 
20 between the ribs and turn sharp laterally as in the holotype (Supplementary Fig. 2 M- 
21 ~~ P). The caudal ribs gradually become short posteriorly and disappear at the 12th caudal 
22 vertebra. The lateral elements of the gastralia are straight, and the median elements are 


23 angulated, each with a single lateral process (Fig. 1 G, H). 


24 The clavicle is relatively broader than that of most pachypleurosaur-like 


13 and nothosaurids", and its expanded posterodorsal portion does not 


25 sauropterygians 
26 have an anterolateral process (Fig. 1 G, H). The clavicle anteriorly contacts the scapula 
27 ~—~*(Fig. 1 G, H; Supplementary Fig. 2 A, B). The scapula develops a distinct ridge along 
28 the midline of the ventral surface and a pit at the posterior end of this ridge. The 
29 coracoid is oval with a coracoid foramen open at the suture with the scapular, 


8.10 nlacodonts?, and even hupehsuchians!> and 
’ Ps 


30 resembling that in saurosphargids 
31  thalattosaurs, but distinctly different from that in any eosauropterygian taxon. The ilium 
32 is not entirely exposed but it is evidently small, much smaller than other girdle elements 
33 (Fig. 1 I, J). In both the new specimen and the holotype, the pubis is plate-like and 
34 extremely round (Fig. 1 I, J; Supplementary Fig. 2 M-P) being reminiscent of that in 
35 saurosphargids and placodonts but different from pachypleurosaur-like sauropterygians 
36 and nothosaurids. An open, drop-shaped obturator foramen is present in the pubis near 
37 its suture with the ischium as in the holotype (Fig. 1 I, J; Supplementary Fig. 2 M-P). 


38 The ischium is kidney-shaped, with an incurved lateral edge and a convex medial 
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1 margin. 
2 The humerus is nearly straight, with a slightly curved posteromedial margin (Fig. | G, 
3H). It is almost uniform, column-like in outline. The ectepicondylar groove is absent, 
4 and an entepicondylar foramen is present. The radius and the ulna are equal in length 
5 and similar in width. The proximal end of the ulna is slightly expanded. There are at 
6 least four roundish carpal ossifications including two large elements, the ulnare, and 
7 the intermedium. The femur is also straight and its posteromedial margin is concave; it 
8 is as long as but its shaft is thinner than that of the humerus (Fig. | I, J). The internal 
9 trochanter of the femur is weakly developed but evident, and the intertrochanteric fossa 
10 is distinct but reduced. Both the tibia and fibula are straight and nearly equal in length, 
11 whereas the former is slightly wider than the latter. At least three tarsals are present, 
12. which are identified as the astragalus, the calcaneum, and a distal tarsal abutting to the 
13 metatarsal I. Both the manus and the pes are small relative to the body size and do not 
14 show a fin-like or a paddle-shaped modification. 


15 Phylogeny. 


16 ‘To investigate the relationship among early sauropterygians, we constructed a new 
17. morphological character matrix, which is largely improved compared to previous 
18  studies!®!* (see Materials and Methods, and Supplementary Information), involving 62 
19 (more than 40 sauropterygian taxa) and 181 characters mostly observed by us on 
20 specific specimens. To further examine whether our new specimen can be referred to 
21. Hanosaurus hupehensis, we coded the holotype specimen of Hanosaurus hupehensis 
22  (IVPP V 3231) and the new specimen (IVPP V 15911) respectively. Results from both 
23 the original matrix of Neenan et al.'°!? and our new matrix indicate that these two 
24 specimens stably form a single clade (Supplementary Fig. 3, 5). Considering their 
25 diagnostic morphology and same locality, we referred this new specimen to the same 
26 species and combined their scores into a single taxon as Hanosaurus hupehensis to 
27 maximize the morphologic information in other analyses here. Using the original matrix 
28 from Neenan et al.'*!°, the basal position of Hanosaurus is recovered (Supplementary 
29 ~+Fig. 4). Using the new matrix, our phylogenetic analyses conduct both maximum 
30 parsimony and Bayesian methods when the latter is employed to Sauropterygia for the 
31 first time (Materials and Methods). The topology generated by the consensus of the 
32 most parsimony trees (Supplementary Fig. 6-8) is fairly comparable to the Bayesian 
33 results (Supplementary Fig. 9, 10). Our phylogenetic results further confirm a 
34 monophyletic clade including all the eosauropterygians, placodontiforms, 
35 saurosphargids, and relevant Triassic marine reptiles (e.g., Helveticosaurus and 
36 Atopodentatus). Here we name this new and previously unnamed monophyletic clade 


37 as Sauropterygiformes (Fig. 2) and Hanosaurus is stably recognized as the basal-most 
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1 taxon of this clade in all most-parsimony trees and Bayesian results. 


4 DISCUSSION 


5 Hanosaurus as the basal-most sauropterygiform. 


6 Considering all the known sauropterygiform taxa on the genus level based on 

7 representative species for the first time, our work illustrates a more explicit phylogeny 

8 of sauropterygiform marine reptiles (Fig. 2, Supplementary Fig. 11). According to our 

9 analyses (see SI for more synapomorphies of main clades), Hanosaurus is recognized 
10  asasauropterygiform reptile based on the possession of: an akinetic palate, a pectoral 
11 fenestration, gastralia segments with five elements, and the crown of marginal teeth 
12 with the concave lingual surface. Moreover, Hanosaurus 1s suggested to be the basal- 
13 most taxon of Sauropterygiformes in terms of the absence of: elongate transverse 
14 _ processes of dorsal vertebrae, radius longer than ulna, and curved humerus that all are 
15 shared by relatively less-basal sauropterygiforms. Additionally, Hanosaurus exhibits 
16 mosaic anatomical morphologies resembling saurosphargids, basal placodontiforms, 
17 _ and basal eosauropterygians. The general body size and shape, the cranial configuration, 
18 and the pachyostosis of the vertebral column, collectively reflect its close affinity to 
19 primitive eosauropterygians including the pachypleurosaur-like forms. However, 
20 compared to any of the known eosauropterygians, the neck of Hanosaurus is shorter 
21 and possesses fewer cervical vertebrae. Given the phylogenetic position of Hanosaurus 
22 asan early-diverging clade, it suggests that a neck of modest length and a cervical count 
23. of about 15 are primitive for Sauropterygiformes, and the neck became even shorter in 
24  Saurosphargidae and Placodontiformes when its length largely increased in 
25  Eosauropterygia, especially in the subgroup, Plesiosauria. Significantly, the disc-like 
26 _coracoid and pubis, and less constricted ischium in Hanosaurus are reminiscent of 
27. saurosphargids, basal placodontiforms, and even earlier deviated thalattosaurians”’ and 
28  hupehsuchians!*. Nevertheless, differing from saurosphargids and placodontiforms, 
29  Hanosaurus also remains a primitive condition, lacking dermal armor, uncinate process 
30 on ribs, and further specialized teeth. The basal status of Hanosaurus and the 
31 phylogenetic result correspond to the fact that it represents one of the oldest members 


32 of Sauropterygiformes that occurred in the Early Triassic”!. 
33 Body plan of early sauropterygiforms 


34 Although as the basal-most sauropterygiform, Hanosaurus does not represent the 


35 terrestrial ancestor of sauropterygians. Apart from being preserved in marine sediments, 
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1 Hanosaurus is certainly an aquatic animal based on its pachyostotic vertebral column 
2 achieving neutral buoyancy in the water, and the reduced ilium with bar-shaped sacral 
3 ribs that could not support the bodyweight to walk on land”. Besides these 
4 modifications, the slender and elongate trunk with the remarkably shortened limbs of 
5  Hanosaurus also featured aquatic adaptations as seen in other marine reptiles**. We 
6 illustrate the body shape of marine reptiles by direct measurements and ratios of various 
7 body regions (see Materials and Methods, and Supplementary Information). Terrestrial 
8 reptiles (e.g., Araeoscelis and Vadasaurus) normally have modest trunk lengths with 
9 considerable sizes of limbs (Fig. 3) for support and movement on land. Conversely, at 
10 the early stage of aquatic adaption in many marine reptilian groups (e.g., Chaohusaurus, 
11 an ichthypterygian; Nanchangosaurus, a hupehsuchian; Pleurosaurus, a marine 
12. rhynchocephalian), a body form combining an elongate trunk and shortened limbs is 
13 consistently present (Fig. 3). The elongate trunk can be additionally reflected by the 
14 higher number of dorsal vertebrae (Supplementary Fig. 12). Many marine reptiles also 
15 have a very long tail, but this long tail seems an ancestral feature inherited from their 
16 terrestrial ancestors” like the early diapsids (Supplementary Fig. 12). 
17. ‘In multiple lineages of marine reptiles, this ‘long-trunk’ convergence reflects the 
18 evolutionary constraints inheriting the axial undulatory movement”? from. their 
19 terrestrial ancestors, and many basal members are anguilliform swimmers such as the 
20 primitive ichthyosauriforms”°*’. Sauropterygia used to be considered as an exception”? 
21 because they were known as paraxial swimmers having relatively short and rigid 
22 trunks”® and using the limbs as thrusters, in which nothosaurs”’ and plesiosaurs are 
23. _—-representatives. However, Hanosaurus possesses more dorsal vertebrae 
24 (Supplementary Fig. 12) and evidently shortened fore- and hind-limbs compared to 
25 other sauropterygiforms, but more similar to other aquatic reptiles like Nachangosaurus 
26 and Chaohusaurus in this aspect (Fig. 3, Supplementary Fig. 12, 13). With the 
27 discovery of Hanosaurus, we herein propose that the basal sauropterygiform marine 
28 reptiles are also dominated by the evolutionary constraint to evolve an elongate trunk 
29 and short limbs at their early stage of secondarily aquatic adaptation. This ‘long-trunk’ 
30 convergent in marine reptiles (Fig. 3, Supplementary Fig. 14) is probably selected by 
31 the coastal marine environment full of rocks and reefs, and/or caused by the mutation 
32 of dorsal vertebra duplication*® that is possibly easier to occur compared to the podial 
33 modifications in limbs. Although the pachyostotic dorsal ribs and the dense gastral 
34 basket preclude flexion of the trunk, the marine reptiles developing pachyostosis could 
35 be anguilliform swimmers*!. We speculate that Hanosaurus and early sauropterygiform 
36 _ reptiles have axial undulatory locomotion at least more than most of sauropterygians if 
37 not as much as that in ichthyosauriforms. Therefore, before evolving the underwater- 
38 flying stroke of pistosaurs, the elongate neck in eosauropterygians, and the turtle-like 
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1 appearance in placodonts and saurosphargids (Fig. 2), Hanosaurus shows an 
2 unexpected ‘long-trunk’ body plan as a possible anguilliform swimmer at the beginning 


of sauropterygiform evolution. 
4 Adaptive radiation of early sauropterygiforms. 


5 After the ‘long-trunk’ convergent evolution of the basal-most members at the very early 

6 stage, sauropterygiforms show a divergent evolution to generate Saurosphargidae with 

7 laterally expanded flat trunks, Placodontia with a body covered with heavy armours 

8 superficially resembling turtles, and Eosauropterygia including abundant 

9  pachypleurosaur-like small members and large plesiosaurs developing long necks and 
10 broad fins (Fig. 2). These different sauropterygiform clades evolved manifold 
11 phenotypes occupying a large region in the morphospace (Fig. 4, Supplementary Fig. 
12 13). In comparison, ichthyopterygians have less change from anguilliform 
13. Chaohusaurus to dolphin-like derived ichthyosaurs, and thalattosaurians show little 
14 bauplan modification besides the changes of head shape and neck length (Fig. 4). Based 
15 on our latest dataset and phylogenetic results of Sauropterygiformes, our tip-dating 
16 result indicates that Sauropterygiform marine reptiles radiated rapidly (Supplementary 
17 _ Fig. 11) from the latest Permian to the earliest Triassic (Fig. 2). High rates are estimated 
18 at the beginning of sauropterygiform evolution, and the three main subclades 
19 (Saurosphargidae, Placodontiformes, and Eosauropterygia) arose within less than 10Ma 
20 from the latest Permian. When more taxa of early sauropterygiforms were added, the 
21 high evolution rates among placodonts*” are not recovered. Another episode of rapid 
22 evolution is at the early stage of eusauropterygian clades leading to the large-sized 
23 nothosaurs and pistosaurs (Supplementary Fig. 11). Apart from the taxon diversity with 
24 more than 40 genera discovered, under our updated phylogenetic framework, 
25 sauropterygiform reptiles took full advantage of the newly emerged marine ecological 
26 opportunities in its early evolutionary burst by evolving novel morphologies with 
27 associated locomotion*’, habitats, and diets. Sauropterygiforms soon rose from the 
28 aftermath of the end-Permian mass extinction to become a predominant reptilian group 


29 conquering the sea. 


31 STAR METHODS 


32 Fossil materials 


33. The new specimen (IVPP V 15911) and the holotype (IVPP V 3231) of Hanosaurus 
34 hupehensis are permanently deposited at the Institute of Vertebrate Paleontology and 
35 Paleoanthropology (IVPP), Beijing, China. As for other observed specimens, we had 
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examined the specimens as many as possible, especially in the richest collections of 
Triassic marine reptiles in China and European countries. We (at least one of the authors) 
examined the true specimens of more than 80% of the total sampled taxa 
(Supplementary Table. 1). The leading author (WW) observed more than 70%. One or 
more of the authors conducted a detailed study on about 50% of the total taxa. To avoid 
bias due to the observation and interpretations of different authors, the state codes for 
all characters were carefully scored by the leading author (WW) with a consultation 


with other authors. As for the genera including several species, if possible, we chose a 


Oo AN DO FP WY FP 


species that is represented by complete, well-preserved adult specimens 


= 
(=) 


(Supplementary Table. 1). 
11 Character matrix construction 


12 ‘To investigate the phylogenetic relationships of Hanosaurus among marine reptiles, we 
13 constructed an updated morphological dataset (62 taxa, 181 characters) focusing on 
14 _ Triassic sauropterygians and relevant reptiles in terms of a series of previous character 
15 matrix and character state coding were largely based on our firsthand observation of 
16 specimens. Compared to the previous most comprehensive dataset of 


16.17 we added 26 new taxa to include all the known placodontiforms, 


17‘ sauropterygians 
18  eosauropterygians, and saurosphargids as the ingroups or the target groups herein, and 
19 sampled the taxa more adequately and representatively than any of the previous 
20 ~ datasets!* to avoid involving redundant species of one genus**. Youngina capensis that 
21 could represent the characters of primitive diapsids was used to root the cladogram. 
22 According to the latest studies from multiple independent research groups!©!7*23>6 
23. and our phylogenetic results using the previous matrix (Supplementary Fig. 3, 4), 
24 — Ichthyopterygia and Thalattosauria have been frequently designated as two successive 
25 sister groups of Sauropterygia. Therefore, with more primitive members of 
26 Ichthyosauriformes**’ discovered from eastern Tethyan fauna, instead of using the 
27. ~whole group of Ichthyopterygia and Thalattosauria as operational taxonomic units in 


28 the previous studies!®!7.!° 


, we used specific species of the basal-most and well- 
29 documented members of ichthyosauriforms (hupehsuchians included) and 
30 _ thalattosaurians as the outgroups to investigate the relationship among sauropterygian 
31 and saurosphargid taxa. Although the monophyly of the marine reptile group combining 
32 Ichthyopterygia, Thalattosauria, and Sauropterygia is controversial, it is beyond our 
33 scope here focusing on the inner relationship within sauropterygians and_ their 
34 relationship with unsolved taxa such as saurosphargids and Atopodentatus. A larger 
35 matrix involving the basal members of Archosauromorpha and Lepidosauromorpha 


35,36,38,39 are better 


36 _ will be considered in our future work when their inner relationships 
37 resolved respectively and more specimens are personally checked by us. The taxa of 


38 stem turtles (or pantestudines) were possibly related to sauropterygians in a few recent 
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1 studies*®*!, and we first-hand observed some specimens‘! and tentatively added them 
2 in our matrix to double test. We added 40 new characters according to the homology 
3 principle and excluded the features dominantly affected by their aquatic lifestyle. We 
4 scored the character states based on the best-preserved adult specimens to decrease the 
5 ontogenetic, sexual, and individual variations. Under these criteria, our character matrix 
6 can be considered on the genus level, the species level, or even the specimen level, 
7 being available when extra operational taxonomic units on any level are added in a 
8 future study. 
9 Phylogenetic analyses and tip dating 
10 To obtain the phylogeny of sauropterygian and allied marine reptiles such as 
11 saurosphargids, we conducted the maximum parsimony inference using TNT version 
12 1.58, All the characters were equally weighted and unordered. The general RAM in the 
13 memory was adjusted to 500 MBytes. Because the number of taxa was fewer than 100 
14 and new technology search was unnecessary, herein we conducted the traditional search 
15 with random additional sequence and tree bisection-reconnection (TBR). The heuristic 
16 search ran 1000 replications of Wagner trees with one random seed and 10 trees saved 
17. per run. We conducted six rounds of maximum parsimony analyses. In the first round 
18 of analyses, we employed the data matrix of Neenan and colleagues!®!”"°, and coded 
19 the character states of the holotype (IVPP V 3231) and the new specimen (IVPP V 
20 15911) separately, when more informative scores from the skull were obtained in the 
21 holotype but more from the postcranial skeleton were coded in the new specimen 
22 respectively (Supplementary Data 1, used with the operational taxonomic unit (OTU) 
23  Hanosaurus hupehensis excluded). In the second round, considering the result 
24 (Supplementary Fig. 3) when two specimens were coded separately in the previous data 
25 matrix, we combined the information into one OTU as Hanosaurus hupehensis 
26 (Supplementary Data 1, used when V 3231 and V 15911 excluded; Supplementary Fig. 
27 ~~‘ 4). In the third round, we coded IVPP V 3231 and IVPP V 15911 separately again in 
28 our updated data matrix (Supplementary Data 2, used with the OTU Hanosaurus 
29  hupehensis excluded) to further test the affiliation of these two specimens 
30 (Supplementary Fig., 5). In the fourth round, the scores were combined as one OTU 
31 Hanosaurus hupehensis (Supplementary Data 2, used when V 3231 and V 15911 
32 excluded) with all 62 taxa involved (Supplementary Fig. 6). Based on the result of the 
33 fourth round, we conducted the reduced consensus tree analysis to find the wildcard 
34 OTUs, which came to be the stem pantestudines and hupehsuchians. Therefore, we 
35 excluded the stem pantestudines in our fifth round (Supplementary Fig. 7), and removed 
36 both pantestudines and hupehsuchians in the sixth round of analyses (Supplementary 
37s‘ Fig. 8). 
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1 To further examine the phylogeny of sauropterygians and saurosphargids, we employed 
2 Bayesian phylogenetic inference using MrBayes version 3.2.7**. Sauropterygians were 
3 all extinct without extant descendants, and our character matrix included only 
4 morphological data, therefore, we used the Mkv model* with variable ascertainment 
5 bias, equal state frequencies, and gamma rate variation across characters and the prior 
6 for the gamma shape parameter was exponential (1.0). For the Markov chain Monte 
7 Carlo (MCMC), two independent runs and four chains, including a cold chain and three 
8 hot chains, were used per run for 10 million generations and sampled every 1000 
9 generations. The average standard deviation of split frequencies (ASDF) was less than 
10 0.01, and the effective sample size (ESS) was larger than 100 both indicating an 
11 adequate result. We conducted two rounds of Bayesian analyses when Hanosaurus 
12. hupehensis was scored based on both two specimens. All 62 OTUs (Supplementary 
13 Data 3) were involved in the first round (Supplementary Fig. 9), and 54 OTUs with 
14 pantestudines and hupehsuchians removed in the second round (Supplementary Fig. 
15 10). 
16 To investigate the evolutionary rates and the divergence times among. all 
17‘ sauropterygiform clades, we performed a Bayesian tip-dating (Supplementary Fig. 11) 
18 using the new data matrix with an unpartitioned model and the constrained parsimony 
19 tree topology in Mrbayes.3.2.7*. We used an independent lognormal“ relaxed clock 
20 model and the Mkv model* with a gamma rate variation*’ across all characters for 
21 likelihood calculation. The parameters of the MCMC process were the same as the 
22 above undated Bayesian phylogenetic analyses. The time tree was assigned a uniform 
23 prior*®, and the root age had an offset exponential prior with a mean of 256.3 Ma anda 
24 minimum of 251.9 Ma. The ages of all marine reptile genera/species in the matrix were 
25 assigned uniform priors respectively with their lower and upper bounds of the 
26 corresponding geologic time stages*” and stratigraphic ranges (Supplementary Table. 
27 ~—-2), except for Psephoderma alpinum being fixed at its last appearance date (201.4 Ma) 
28 because of having the youngest last appearance date among all taxa we considered. 
29 Specimen measurements and body shape analyses 
30 We assembled a series of linear measurements from representative marine reptiles 
31 (Supplementary Data. 4). Our dataset includes all sauropterygiform subgroups and 
32 other major lineages of Mesozoic marine reptiles such as Hupehsuchia, Ichthyopterygia, 
33 Thalattosauria, marine lepidosaurs including marine sphenodontians and mosasaurs. 
34 We sampled and measured 40 nearly complete skeletons of 40 species/genera 
35 representing basal to crown taxa of each group (Supplementary Data. 4). According to 
36 our phylogenetic results and literature, we reconstructed a cladogram of all these 
37 sampled marine reptilian taxa (Supplementary Fig. 12). Measurements of the dataset 
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1 were made either using a vernier caliper and flexible rule directly from the specimens 
2 or using the ImageJ program from the published photos and corresponding scale bars 
3 in literature. Our dataset included nine discrete and unrelated linear values 
4 (Supplementary Data. 4): the skull length (SL) from the anterior margin of the snout to 
5 the posterior end of the occiput, the skull width (SW) as the maximum width of the 
6 dermal cranium, the neck length (NL), the body length (BL) (same as the trunk length 
7 inthis study) from the first dorsal to the last sacral vertebra, the body width (BW) (same 
8 as the trunk width in this study) measured from the maximum width of the rib cage 
9 when the specimen was preserved in dorsal view or the length of the longest rib in 
10 lateral view, the tail length (TL), the tail width (TW) as the maximum width of the 
11 proximal region of the tail in dorsal view or the width of the largest caudal centrum 
12 including the length of the associated transverse process, the forelimb length (FL), and 
13 the hind limb length (HL). 
14 Our analyses were conducted using the paleontological statistics software package 
15 (PAST) version 4.0°°. To test the convergence of long trunk and short limbs, we first 
16 conjunctively considered these two features and made a binary plot using the 
17‘ dimensionless ratios of BL/FL and BL/HL as the coordinates with all original 
18 measurements (BL, Fl, and HL) already logio-transformed (Fig. 3). Secondly, we made 
19 a bar graph based on the ratio of BL/BW to show the elongation of the trunk itself 
20 (Supplementary Fig. 12). Thirdly, we created a stacked bar chart of NL/BL/TL to 
21 compare the percentages of these three parts in different taxa (Supplementary Fig. 12). 
22 These results should be comprehensively considered, because some features might be 
23 ignored due to others, such as that the BL of Henodus (a highly specialized placodont) 
24 occupied a large percentage in the entire length but its trunk was not actually elongated, 
25 while the body of Pleurosaurus (a marine sphenodontian) was extremely long but not 
26 that notable with a tail even longer. These nine measurements could outline the body 
27 _ shapes, therefore, we collectively used aforementioned nine linear measurements being 
28  logio-transformed and employed principal component analysis (PCA) to demonstrate 
29 the shape variances of the sampled marine reptiles. We also conducted principal 
30 coordinate analysis (PCoA), which showed almost the same results as PCA herein. The 
31 principal component | axis (PC1) reflected about 86% of the variance showing positive 
32 eigenvector coefficients with all the nine measurements, hence PC1 seemed mainly 
33 related to the size with large marine reptiles on one side and small ones on the other 
34 (Supplementary Fig. 13). PC2 described more than 7% of the variance being evidently 
35 negative (-81%) with the neck length, which explained the elongation or shorten of the 
36 cervical region with Plesiosaurus and Cartorhynchus (an ichthyosauriform) at the two 
37 opposite ends. PC3 explained over 3% of the variance, and its eigenvector coefficients 
38 were mostly negative (-64%) with the body width and mostly positive (61%) with the 
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tail length when Henodus and Anshunsaurus (a thalattosaurian) were at the two 


bh 


extremes respectively. We preferred the binary morphospace using PC2 and PC3 to 


N 


3. demonstrate the shape disparities of multiple marine reptilian lineages (Fig. 4). 
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Figure Legends 


Figure 1. The first complete skeleton of Hanosaurus hupehensis. VPP V 15911. 


(a-b), new specimen preserved in ventral view; (c-d), skull; (e-f), the medium region of 
cervical vertebrae; (g-h), the right side of the anterior region of the trunk with the 
pectoral girdle and the right forelimb, the line drawing is a mirror of the opposite side 
along the midline; (1-j), the left side of the sacral region and the anterior region of the 
tail, with the pelvic girdle and the left hind limb, the line drawing is a mirror of the 
opposite side along the midline. (scale bar=5 cm in a-b, =1 cm in all the others; see 


abbreviations in Supplementary Information) 


Figure 2. Dated cladogram of Sauropterygiformes based on the analyses of our 


new dataset on Triassic marine reptiles. 


Taxonomic groups are marked by different colors. Topologies are combined with 
results from both maximum parsimony and model likelihood (Supplementary Fig. 5- 
10). Colored bars of each taxon indicate the duration of discovered fossil remains. 
The ages of the nodes are obtained from the posterior median evaluations 
(Supplementary Fig. 11). Geologic time follows the latest monographs (Materials and 
Methods). 


Figure 3. Convergence of body plan in secondarily aquatic reptiles reflected by 
the ratios of trunk/forelimb length and trunk/hind limb length. 


Taxonomic groups are labeled in different colors, and relative phylogenetic stages 
(basal, medium, derived) are plotted using varied shapes. Terrestrial, coastal marine, 
and pelagic marine taxa are respectively clustered in distinct areas. From basal to 
derived members, Sauropterygiforms and ichthyopterygians show a trend from axial 
to paraxial locomotion, while marine sphenodontians and mosasaurs are from paraxial 
to axial. This diagram illustrates the body plan convergence of long trunk and short 


limbs at the early aquatic adaptation stage. 
Figure 4. Morphospace based on principal component (PC)2 and PC3. 


Body shape principal component analysis reflects the morphologic variances among 
different Mesozoic marine reptiles. Black lines outline the area occupied by 


sauropterygiform reptiles. 
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